Electromechanical Behavior of BaTiOa from First Principles 
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Using an effective Hamiltonian parametrized from first-principles, Monte Carlo simulations are per- 
formed in order to study the piezoelectric response of BaTiOs in the ferroelectric tetragonal phase as 
a function of temperature. The effect of an electric field on the phase behavior is also illustrated by 
a simulation of the transformation of a rhombohedral domain into a tetragonal one under a strong 
field. 
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One of the earliest technological applications of ferro- 
electric materials was in the area of electromechanical 
transducers, the physical basis of which is the piezoelec- 
tric effect. The piezoelectric coefficients characterize the 
(linear) change in polarization in the presence of an ex- 
ternal stress, or equivalently, a change in shape under 
the application of an external electric field. While a 
phenomenological, descriptive framework for such elec- 
tromechanical effects has been available for a long time, 
a quantitative microscopic understanding of the response 
of individual materials is lacking, and further technical 
development has to proceed largely by trial and error. 
The importance of obtaining a deeper understanding has 
been highlighted by recent reports [Q of giant piezoelec- 
tric response in single crystals of relaxor perovskites of 
the form PMN-PT and PZN-PT. The structural com- 
plexity of these materials and the variety of phenomena 
that might be involved make it rather difficult to offer a 
clear understanding of their observed properties. 

It is in this context that first-principles calculations 
can help. They can be used to selectively "turn off" 
features of the system and study its response in condi- 
tions that are very difficult or impossible to realize in 
the laboratory. Besides, they can provide a microscopic 
view of the materials properties. So far, the primary uses 
of first-principles calculations relevant to this area have 
been calculations of piezoelectric coefficients for simple 
materials at zero temperature [0. In this work, as a step 
towards the theoretical treatment of the electromechan- 
ical response of complex systems, we present the first 
calculations of the piezoelectric response as a function of 
temperature for BaTiOa, an important member of the 
perovskite family of ferroelectric compounds. We also 
illustrate the infiuence of electric fields on the phase be- 
havior of BaTiOa. We base our approach on a scheme 
which has proven very successful yj in the description of 
the rich phase diagrams of perovskite oxides: an effec- 
tive Hamiltonian which contains the physically relevant 
degrees of freedom of the structure (notably the "soft 
mode" ) is constructed on the basis of high-quality first- 



principles calculations, and the statistical mechanics of 
the system is then studied by Monte Carlo simulation. 
The usefulness of the method is not restricted to the 
calculation of piezoelectric coefficients and the study of 
generalized phase diagrams. Indeed, the elastic and di- 
electric response can be similarly computed, and it can 
form the basis for an analysis of non-linear effects Q . 

The appropriate thermodynamic identity for a crystal 
in the presence of a (possibly anisotropic) stress and an 
electric field is [pi 



dU = TdS + a^d-q^ + E,dD, 



(1) 



Here U is the internal energy of the crystal (per unit 
volume), T and S are the temperature and entropy, Ui, 
and rji, are the components of the stress and strain tensors 
in the Voigt notation , i?i is the jth component of the 
macroscopic electric field, and Di is the corresponding 
component of the electrical displacement vector (in SI 
units, D = eoE-l-P). 

We parametrize the energy U of the system by means 
of an effective Hamiltonian iiTeff which represents a Tay- 
lor expansion of the energy surface around the high- 
symmetry cubic perovskite structure. i?eff is written in 
terms of the dynamical variables which are relevant to 
the low-energy distortions: the amplitudes {u} of the lo- 
cal modes (three degrees of freedom per unit cell) which 
represent the "soft" transverse optical phonon and are 
directly related to the polarization of the crystal [P = 
{Z* /V) ^ u, where Z* is the mode effective charge and 
V is the cell volume] ; a set {v} of displacement variables 
representing the acoustic modes; and the six components 
of the homogeneous strain ry. The parameters of the 
energy expansion, including the on-site local-mode self- 
energy, the interaction between local modes (both short- 
range and dipole-dipole), the elastic energy, and the local 
mode-elastic coupling, are computed using highly accu- 
rate first-principles LDA calculations with Vanderbilt ul- 
trasoft pseudopotentials . More details about the con- 
struction of the effective Hamiltonian are given in Ref. [Q, 



where the method is shown to provide a good account 
of the phase transition sequence in BaTiOa within the 
hmitations of an approach based on low-energy distor- 
tions. The extension of the standard MetropoUs Monte 
Carlo algorithm to include the effects of stress and elec- 
tric field involves replacing the Boltzmann probability 
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factor exp(-/3;7^) by exp[-/?([/^ 
acceptance criterion for state j . For a given temperature, 
stress, and field, the strain 77 and the mode variables are 
allowed to fluctuate, their average values determining the 
shape and net polarization of the system. 

As an important application of the method, we com- 
pute the piezoelectric response of the tetragonal (ferro- 
electric) phase of BaTiOs (point group Amm), which is 
stable from approximately 278K to 403K and exhibits a 
spontaneous polarization that we take to be along the z 
axis. The relevant coefficients are 



dw = 



dE,, 



da. 



/3(AF,A77,) 



(2) 
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where /S.X — X — {X) and the averages are computed 
using the extended Boltzmann factor defined above [gj. 
These equations suggest three different ways to calculate 
the response: direct calculations of the average strain 
as a function of applied field, or the average polariza- 
tion as a function of (anisotropic) stress, and compu- 
tation of the statistical correlation between polarization 
and strain. The latter is conceptually the simplest, al- 
though relatively long simulations (on the order of 100000 
Monte Carlo sweeps (MCS) 0]) are needed to obtain 
good statistics. Of the direct approaches, the field-strain 
calculations are the most efficient as only one series of 
calculations for varying E^ is needed to compute the 
most common coefficients d^i and ^33 (representing re- 
spectively the transverse contraction and the longitudinal 
expansion under the application of a field parallel to the 
ferroelectric axis). 

An effective Hamiltonian based on a finite Taylor ex- 
pansion of the energy should not be expected to repro- 
duce perfectly the behavior of the material at relatively 
high temperatures. In particular (see Ref. Q]) the theo- 
retical transition temperatures are progressively shifted 
downwards with respect to the true ones p4|| . In order to 
provide a better comparison of our results for ^31 and ^33 
to experiment (Fig. n^), we have therefore rescaled linearly 
the theoretical temperatures so that the end points of the 
range of stability of the tetragonal phase coincide. The 
agreement with the available experimental data is very 
good, and the general trend of the temperature depen- 
dence (upward tendency in ^33, downward in ^31) corre- 
sponds to the pseudo-divergent behavior observed exper- 
imentally near the tetragonal-to-cubic transition |l3] . In 
the thermodynamic limit, the correlation and field-strain 
approaches should be completely equivalent. Since we 
use finite runs and finite simulation boxes (12 x 12 x 12 or 
1728 unit cells for the correlation analysis and 10 x 10 x 10 
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FIG. 1. Piezoelectric coefficients as a function of tempera- 
ture in the tetragonal phase. Open triangles and circles repre- 
sent values computed from the strain response to an electric 
field and from polarization-strain correlations, respectively. 
Sohd squares show experimental data points (Ref. |J10|). 

or 1000 unit cells for the field-strain calculations), the re- 
sults are similar but not identical. 

Our method can also be applied to the study of the 
stability of the different phases as a function of the ex- 
ternal field and stress. As an illustration of the influence 
of electric fields on the phase diagram, we performed a 
simulation of the effect of a strong field on the strain and 
electrical polarization of a BaTi03 sample in the rhom- 
bohedral (R) phase. This phase, with point group 3to, is 
stable at temperatures below 183K, and exhibits a spon- 
taneous polarization along one of the original (111) cubic 
directions, and a corresponding strain deformation with 
respect to the cubic phase satisfying rji = 772 = 773 and 
^4 = % = %■ Choosing a simulation temperature of 
lOOK, we applied an electric field along the cubic z axis, 
progressively increasing the strength of the field up to 
a maximum value of 5 x 10'^ kV/cm. The evolution of 
the strain and polarization as a function of the field is 
shown in Figs. and pi respectively. The rhombohedral 
symmetry is immediately broken by the field, as can be 
seen by the iji = r/2 ^ rjs splitting (and the analogous 
one among the polarization components). In the approx- 
imate range 2 — 2.5 x 10'^ kV/cm, there is a noticeable 
split in 771 and 772 , all the off-diagonal components of the 
strain save 7/4 (27723 in the standard tensor notation) fall 
to zero, and Pi falls to nearly zero. Thus, this range cor- 
responds to an orthorhombic phase with the polarization 
in the yz plane (different simulations "pick" at random 
between the yz and the xz planes). For larger fields, 
the equality between 771 and 772 is restored, all the off- 
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FIG. 2. Strain vs. field for the phase transformation from 
rhombohedral to tetragonal under a strong field. Off-diagonal 
strain components are magnified for clarity. 



diagonal strain components fall to zero, and P2 follows 
P\ in its earlier drop. The resulting phase is tetragonal, 
with P along the [001] direction. The behavior emerging 
from our simulation parallels that seen experimentally |l|] 
in the PMN-PT and PZN-PT giant-response materials: 
individual domains are in a rhombohedral structure, and 
as an electric field is applied and progressively increased 
along z, the polarization rotates from a [111] direction to 
a [001] direction. 

In conclusion, we have shown how the electromechan- 
ical response of a system can be computed from first 
principles using an effective Hamiltonian suitably aug- 
mented by terms representing the influence of applied 
(anisotropic) stresses and electric fields. As an appli- 
cation, we have presented the first calculations of the 
piezoelectric response of ferroelectric tetragonal BaTiOs 
as a function of temperature. We have also performed 
simulations illustrating the effect of electric fields on the 
phase behavior of this material. 
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FIG. 3. Polarization vs. field for the phase transformation 
from rhombohedral to tetragonal under a strong field. 
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